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1. The Gribov no-pole condition 

In the seminal work [|l|] Gribov noticed that neither the (covariant) Landau nor the (noncovari- 
ant) Coulomb gauge are suitable to select a single gauge-field representant on each orbit of gauge- 
equivalent fields. Later on, it was recognized that this so-called Gribov ambiguity is a generic 
problem for non-Abelian gauge fixing [^. 

The main observation by Gribov is that, upon considering an infinitesimal gauge transforma- 
tion parameterized via one observes that the Landau-gauge condition, d^A^ = 0, 
is also fulfilled by the equivalent field, i.e. A'^, d^A'^ = 0, if the Landau-gauge Faddeev-Popov op- 
erator 



has zero modes, since we must have, at first order in CO, 



(1.1) 







d^A^ 



(1.2) 



If one removes the possibility of such zero modes by modifying the path integral measure, at least 
infinitesimally-related gauge copies will be ignored in any computation of expectation values. Let 
us also notice here that the existence of zero modes of the Faddeev-Popov operator invalidates, 
strictly speaking, the mathematical correctness of the textbook Faddeev-Popov quantization. In- 
deed, the latter procedure is based on the appropriate functional generalization of the identity 



5/ 



x=x*J{x*)=y 



Jdx5[fix)-y]g{ 



x=x'' ,f{x*)=y ' 



(1.3) 



where ^ 



plays the role of the Faddeev-Popov determinant. However, the above expression 



df 



is only valid if f{x) = y has a single solution at x = x*, thereby also assuming that 
Additionally, assuming a positive Faddeev-Popov determinant also allows to drop the absolute 
value operation, something which is usually done. If /(x) = would have multiple solutions, let's 
say at X = X* , one should use 



Jdx5[fix)-y]g{x)=l^Lix) 



df 



dx 



(1.4) 



x=x* 



but, to our knowledge, the functional generalization of this expression cannot be brought into a 
workable action for a gauge fixed SU(N) Yang-Mills theory, mainly because (1) a single Faddeev- 
Popov determinant can be lifted up into the action by the introduction of ghost fields, but a sum of 
determinants does not allow for this trick, and (2) the sign of the determinant should be taken into 
due account. 

So, if one desires to stay as close as possible to the conventional Faddeev-Popov procedure, 
one better restricts the gauge-field region, over which the path integral runs, to the following set of 
configurations: 

a = {a';^{x) : d^A'l^ix) =0,^"''{x,y) > o} . (1.5) 
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This Gribov region, with the (first) Gribov horizon being its boundary dO., exists out of the relative 
minima under gauge transformations of the functional (^[A] = / d^xA^. Gribov considered the 
ghost propagator written as 

By imposing 

o{k^)<\ for k^>Q (1.7) 

he required ||l|] that the ghost dressing function k^'i^{k^) cannot have a pole at finite nonzero mo- 
menta. Since the ghost propagator is actually given by 

it is clear that the inequality ([L7|), known as the no-pole condition, should be equivalent to the 
restriction of the functional integration to the Gribov region Q.. 

In this proceeding, based on the work [^], we shall examine, as general as possible, the con- 
sequences of the no-pole condition; in particular, we shall focus on the specific case that the num- 
ber of space-time dimensions is 2. Why our interest in this case? At present day, there is little 
doubt that lattice simulations of the Landau gauge prefer the so-called decoupling/massive solu- 
tion, in which case the gluon propagator attains a finite nonzero value at zero momentum, while 
the ghost propagator strives to a Xjp^ singularity at low momenta; this at least in 3 or 4 space- 
time dimensions [@, H |[ 0, |, |, |l^]. This was also discussed from the analytical viewpoint in 
e.g. Ill, 12, 13|, |l4[ [l5|, l^]. These results replaced the for a while rather popular scaling scenario. 



in which the gluon propagator vanishes at zero momentum, while the ghost propagator becomes 
more singular than 1/p^ [ p^ , p7| ]. The physical relevance of this propagator research program lies 
in their use in phenomenological applications, next to trying to gain a better understanding of con- 
finement, or chiral symmetry breaking, in the Landau gauge. We refer to several other contributions 
in this proceedings to illustrate this. 

In contrast with the J = 3,4 cases, the scaling scenario keeps being supported by the d = 1 
lattice data, as is visible from e.g. [Q, [l^, 19]. A careful analysis based on rigorously derived 



bounds allowed to show that ^(0) = persists in the infinite-volume limit []g] (see also []19|]). 
Here, we present a summary of the detailed analysis presented in [^]. The aim is to unravel a 
theoretical reason behind the impossibility to find a gluon propagator with ^(0) > in (i = 2. 
We shall see that the problem with a massive kind of gluon propagator always reveals itself under 
the form of a logarithmic infrared singularity in the ghost self-energy, viz. the no-pole function 
o(k^). Evidently, such singularity should be excluded, as it clearly violates the Gribov condition^ 
(pT?!). For completeness, we also performed a similar analysis in J = 3,4, in which case we are not 
forced to conclude that i^(0) = 0, meaning that a decoupling/massive gluon is perfectly possible 



in fif = 3,4, and not necessarily in contradiction with the requirement (1.7) 



In any case, we do not expect to find a singular term in the quantum correction a(]c 
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2. Scrutinizing the d = 2 ghost propagator 

2.1 Infrared singularity in and bounds on o(k^) using one-loop perturbation theory 

As a warming up exercise, let us first consider the one-loop-corrected ghost propagator in the 
Landau gauge, which reads 



1 



{k-qf 



(2.1) 



with 5"'' &{q^) ^^v(^) the Landau-gauge gluon propagator and ^^v{q) = {S^v — <iii<lv/q^) the 
transversal projector. For the moment, we also assume a tree-level ghost-gluon vertex igf"'^'^ky. 



The expression (2.1) can be simplified and rewritten into the following form 

1 1 



k^l-a{k'^) 

where CJ(A:^) is the momentum-dependent function given by 

2,T^^^v f d'^<7 1 



(2.2) 



{InY {k-qf 



(2.3) 



This corresponds to the usual resummation of an infinite set of diagrams into the IPI ghost self- 
energy. Notice that this resummation only makes sense when (y{k^) < 1, i.e. when the no-pole 



condition ( |1.7| ) is satisfied. 

We can then manipulate^ by two consecutive partial integrations 



o{k^) _ 

g^N ~ Stt 

1 

~ 871 

_ 1 

~ 871 

~ Stt' 



+ In (x) 

- In (k^) &{k^ 

- In {k^) S){k^ 



D{k^ 


)-m 






D{k^ 


)-D(0) 




k2 


b{k^ 


)-D(0) 



dxln(x) Si'{x 

dxln(x) Si'{x) 
xln(x) -x] Si'{x) 



k^ 



+ J dx [x In (x) 



1 J D{k^)-D{0) 



Btt 



k^ 



In {k^) &{k^) + [k^ In [k^) - k^] &'{k^^ 



+ dx[xln(x) -x] 9"{x) 



^We refer to [HI for an alternative derivation. 
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where we assumed^ for simplicity that &{x) ~ 1/x at large x and that also & {x) goes to zero 
sufficiently fast at large momenta, e.g. as 1 /x^. We have denoted with D{k^) a primitive of &{k^). 
In the limit — > we then obtain 



- \ ^(0) - lim In ik^) &(0) + f dx[xln (x) - x] 2" {x) \ . (2.4) 

iTT I ^2^0 ' JO I 



Clearly, we encounter a small-momentum singularity proportional to — ^(O)ln(^^), which can 
only be avoided upon demanding that ^(0) = 0. This is the first indication that it might be difficult 
to have a decoupling/massive gluon, S>{Q) > 0, in 2 space-time dimensions. In specific cases, 
the infrared sickness oi d = 2 Landau-gauge theory with a massive gluon was already noticed in 
[20, The preference for the scaling solution in <i = 2 was very recently also proclaimed from 



a somewhat different viewpoint in [P4|. 

2.2 Infrared singularity in and bounds on o{k^) using the ghost Dyson-Schwinger equation 

Since the analysis so far is purely based on one-loop perturbation theory for the ghost propa- 
gator, albeit with a completely general gluon propagator, we need to certify our conclusion beyond 



this approximation. Therefore, we shall first derive a set of bounds for the one-loop result ([2.3|), 
valid for general d. We need that (see ||3|, |23| ]) 



+ ^ ^ ^ " ^1^'' 1 + d/2;q^/k^) e{k^ - q^) (2.5) 



where the Gauss hypergeometric function 2Fi{a,b;c;z) is defined for \z\ < 1 by the series 
17 ( h \ f {b)n z" .ab a{a+l)b{b+\) 

with {a)n = r{a + n) /r(a) the Pochhammer symbol. The previous series is converging for c / 
0,-1,-2, . . ., and in addition also for |z| = 1 if ^e{c — a — b) > 0. In the latter case, we have for 

1) = ^!^^%^. (2.7) 

Y{c — b)Y{c — a) 
We also introduced the usual solid angle integral 

f In'^l'^ 
Q.d= dHrf = . (2.8) 

For d = 3,4 and with z G [0, 1], it can be shown that the hypergeometric function 2^1 (1, 1 — <i/2; 1 + 
d/2;z) has its maximum value, equal to 1, at z = 0, and its minimum value, equal to d/ {2{d — 1)), 



■'More mathematical details, as well as a more profound discussion on the assumptions made at the level of the 
large/small momentum behaviour of the gluon propagator can be found in [b||. 
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at z = 1. For d = 2, the function 2^1 (1, 1 —d/2; 1 +d/2;z) = 2^1 (l,0;2;z) is identically 1 for 
|z| < 1. As such, for general d, we can propose the following estimate 



2{d-l) - g^N 



with 



Id 



d-i 



{2nY d Jo 
for the one-loop ghost self-energy 



d{k^-q^) e{q^-k^^ 



+ 



(2.9) 



(2.10) 



1 — COS^(0i 



' cos(0i 



Q.r, d-l 



{2nY d Jo 



dqq''-^ ^(^2 

d{q^-k^ 



d{k^-q^^ 



■2Fi{l,l-d/2;l+d/2;q^/k^) 



+ ■ 



■2fi(l,l-^i?/2;l+J/2;A:V^^^ 



(2.11) 



From the bound ( |2.9[ ) we derive an equality for J = 2, in which case there is a singularity, since, as 
already discussed in the previous section, /j displays one. 

To go beyond perturbation theory, we take a look at the Dyson-Schwinger equation for the 
ghost propagator [11, |l^, Q for the d = 2 Landau gauge 



1 



k^-g'-N J ^p,rUk,q)ky&{q^)P^y{q)^{p^) , (2.12) 



where p = k — q and where the gluon and the ghost propagators, &{k^) and W{k^), are now fully 
dressed propagators. The quantity igf"'^'''k^r^y{k,q) corresponds to the full ghost-gluon vertex. 



Upon using eq. (2^), we can consider 

= §^ / 702PAr,v(^,^)^v^(^^)P,v(g)-^ ^_^. 2) • (2-13) 



For a tree-level ghost-gluon vertex r^y{k,q) = d^y and using the transversality of the gluon prop- 
agator we write 



d^q 

(2 



1 1 

p^ 1 - o{p^) ' 



(2.14) 



or, using a suitable base for polar coordinates, 
o{k^) rqdq 



47r2 



_ 1-COS2(0) 

ad 



(2.15) 



with p2 =k'^ + q'^ -2kqcos{d). 
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Let us first assume there is no ghost enhancement, i.e. that a(0) < 1. In such case, we easily 
derive a lower and an upper bound for the l.h.s. of Eq. ( 2.15 ) as follows"^: 



I < \ < 



g^N - l-a(O) 



where 



' qdq 
47r2 



'^^ l-cos2(0) 

do 



(2.16) 



(2.17) 



This integral / is the one-loop ghost-self-energy contribution we already encountered before, so 
we know that / will display an infrared singularity proportional to — \im^.2^Q\n{k^) md = 2 unless 
^(0) = 0. A fortiori, the o{k^) of eq. ( 2.15 ) has a logarithmic infrared singularity, unless ^(0) = 0. 

A little more care is needed when a(0) = 1. Let us then write \ — o{k^) ck^'^ for small k^ 
and analyze 



' qdq 



+ 



4k^ 
qdq 
An 



&{q^ 



^ 1-COS2(0) 

dd - 



2 



2^ _ 1-COS2(0) 

dd ^ 







\ — o{p^) cp 



2e 



(2.18) 



The quantity in square brackets in the last integral is finite at = if the behavior of o{k^) is given 
by 1 — ck^'^ + ^{k!') with b > 4e. This quantity goes to 1 at large momenta and its absolute value 
is bounded from above by some positive constant M if (y{k^) G [0, 1]. As such, we majorate as 
follows 



< 



' qdq 
4k^ 



&{q^ 



2^^^^ 1-COS2(0) 

dd , ^ +M/, 



cp 



2+2e 



(2.19) 



with / the same integral as before. For the other integral we can use the general result ( |2.5[ ) for 
d = 2 and v = I + e, obtaining 



qdq 



4712 

fl^qdq &{q^) ^ , 2 /,2n ^ T ^ ^ ^ , 2^ 



2^ l-COs2(e) 

de- 



ep 



2+2e 



qdq Sl{q^ 



From the series representation (2.6) it follows that the derivative 2p\ («, b; c;z) can be brought in the 
following form: 

^2Fi{aMc;z) = —2F\{a + \,b+\;c+\;z) . (2.21) 
oz c 



In particular. 



— 2F\{'\-+e,e;2;z) 
az 



[\+e)e 



2Fi{2 + e,\+e;3;z) 



(2.22) 



i.e. it is positive, implying that this hypergeometric function attains its largest value at z = 1. As 
9?(c — a — b) = \— 2e,\Ne. have that for e < 1 /2, the hypergeometric function is finite at z = 1, with 



■^Here, we assume that o{k^) attains it maximum value at = 0, which was already done by Gribov y]. In some 
cases this can be shown (jsjl, but going into those details would lead us too far now. In fact, even if a(A^) < 1 would be 
maximal at =kl > 0, the following argument would still work out by replacing tT(0) with o{kl). 
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a value given by r(l — 2e)/ [r(2 — e)T{\— e)]. Thus, for e < 1/2 a singularity can only appear for 
0. Using the trapezoidal rule, 



' dxf{x) = ^ [/(^) +/(«)] + ff{b-a)' 



we subsequently find in the limit ^0 

k^ &{k^)2Fi{l+e,e;2;l) 



lim 



Stic 



k 4kc q 



(2.23) 



(2.24) 



As before, we note the presence of an infrared singularity, due to the last term, whenever ^(0) > 0. 



Both terms on the r.h.s. of the inequality ( |2.19| ) develop a — ln(0) singularity proportional to ^(0). 
Since the lower bound in eq. ( |2.16| ) always applies to eq. ( 2.15| ) and / displays the same singularity 
as obtained for the upper bound, we are forced to conclude once more that o{k^) can be finite at 
k^=0 only if ^(0) = 0. 

To arrive at the general conclusion that in = 2 we must have ^(0) = in order to avoid 
an infrared singularity in the ghost self energy, we had to assume that o{k^) = \ — ck^" + i?{k^) 
with \>2e and b >4e in the case ghost enhancement occurs. The d = 2 lattice data [ p^ suggest 
for the ghost propagator an infrared exponent e ^0.15, in reasonable agreement with the scaling 
solution of [23, 2^, 26] that predicted e = 1/5. Both lattice and analytical estimates are such that 
the condition 1 > 2g is satisfied. The demand on the subleading exponent, b > 4e, is something we 
must assume, since there are neither lattice nor analytical estimates for this guy. 

The attentive reader shall have noticed that so far, we relied on the use of a bare (tree level) 
ghost-gluon vertex. Such approximation is usually made in Dyson-Schwinger studies. Though, 
to obtain a qualitatively precise comparison between the lattice and analytically-obtained Dyson- 
Schwinger ghost propagator, one should go beyond a tree-level vertex [p7|]. Fortunately, one rec- 



ognizes that, by taking in Eq. (2.13) the full vertex rxy{k,q), instead of the tree-level one S^y, the 
above derived results still apply, as long as the vertex itself remains infrared finite, as seems to be 



confirmed by lattice data [28, 2^. This illustrates the power of the presented bound analysis. 



3. What about J = 3 or J = 4? 

Having done all the efforts to analyze the d = 2 ghost propagator from its Dyson-Schwinger 
equation, one could wonder if not a similar result could be derived in J = 3 or J = 4 and, if so, if 
we can draw any conclusions as in J = 2? For the sake of presentation, let us stick to the d = 4 
case in this proceeding; the d = 3 case is completely analogous. The following analysis is again 
originating from [||]. 



The Dyson-Schwinger equation for o{k ) is given by [g4p 



a{k^ 



1-2^3 + ^ig'N / dq 



{2k) 



■&{q') / da. 



1 — COS^(0i) 



(3.1) 



where i?3 and are the renormalization constants for the ghost propagator and the ghost-gluon 
vertex respectively. We have again opted for a tree-level ghost-gluon vertex here, but as for the 
d = 2 case, if the vertex is infrared finite all results generalize to a nontrivial vertex. 
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Thanks to the the non-renormaUzation of the ghost-gluon vertex in Landau gauge [30|, we 
may set ^\ = \. We shall build the argument on 

■--(^')- r,,|;;l^(,.)/,^, ■ (3.2, 



for d general, in the spirit of dimensional regularization. We will show that no infrared singularities 
occur for 2 < J < 4, consequently the ultraviolet divergences when d ^ A we. treated by the iF- 
factor(s). We again derive bounds, but this time for ad{k^)- We have 

g^N - 1 - a,(0) ' ^'-'^ 

ifOdik^) < OdiO) < 1, and 

</(e)+M/(0), (3.4) 
if 04{k^) < OdiO) = 1 with Gdik^) ^ 1 - ck^" for small k^. We introduced 



The basic inequalities ( [2.9| ) can again be used to write down that 

Odik^) / Id 



iiOdik'^) < GdiO) < 1, and 

^<Iie)+MId, (3.7) 

if Od{k^) < Od{0) = 1. The integral /j has been defined already in eq. ( 2.10| ), making clear that Id 
for 2 < J < 4 is finite, even if ^(0) ^ 0. We can also evaluate /(e) further, yielding 

I[e) = / , 2Fi[^+e,l+e-d/2;l+d/2;q /k ) 



{2kYc T(d/2) k^ 

Ife<<i/2— 1, the hypergeometric function attains it maximum value, 1, at z = 0, while the above 
integrals are all well-behaved for k^ — )■ 0, as we assumed d >2. If, on the other hand, e > d/2 — 1, 
the hypergeometric functions' maximum is at z = 1, with convergence if e < (J — l)/2. We can 
thus state that 

^ < 0.9) 

for some M' > 0. Hence, for <i = 4 we have a finite aj(0) if i^(0) is also finite, but it does not have 
to vanish, in contrast with our d = 2 findings. 

For completeness, we should mention that when e > d/2 — I, in which case the hypergeo- 
metric function 2F1 {I + e, \ + e — d /2; \ + d /2;z) would not be convergent at z = 1, we are unable 
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to derive a sensible upper-bound^. Such values for the exponent e would correspond to a ghost 
propagator which would be very infrared singular, in particular for d = 4, the ghost would need to 
be at least as singular as \/k^ . Such a strong singularity has never been proposed to our knowledge, 
neither numerically nor analytically. 



4. What do our d = 2 results imply? 

We have thus analytically shown, using as sole input the ghost Dyson-Schwinger equation and 
the Gribov no-pole condition, that in <i = 2, one can only find a gluon propagator that vanishes 
at zero momentum, ^(0) = 0. If ^(0) > 0, a logarithmic infrared singularity occurs in o{k^), 
which makes it impossible to fulfill the Gribov no-pole condition. The latter serves as a diagnostic 
tool to check whether one did not cross the first Gribov horizon, where the Faddeev-Popov oper- 
ator vanishes. The used mathematical tools are easy to follow step by step, and do only rely on 
some general properties of e.g. the gluon propagator, rather than on (un)controllable approxima- 
tions. Our results in a way provide the analytical "reason" behind the d = 2 large-volume lattice 
results and their accompanying infinite- volume limit [Q, ^, [I^, |l^]. Indeed, there is no sign of a 
decoupling/massive scenario using lattice simulations of <i = 2 Landau-gauge theory. 

Our results indicate that the finding of the decoupling/massive solution, with ^(0) > 0, as 



reported in [25, 26], should be reconsidered. Moreover, we believe they also pose a challenge 



for the interpretation presented in works like [31, 32, ^] (see also [ p^ for a critical assessment): 
using the fact that one can find an infinite number of solutions to the Dyson-Schwinger equations 



[ ]16| , [34| ], it was (and still is) speculated that choosing one of these solutions would correspond to 
a supplementary, nonperturbative gauge-fixing choice, on top of the underlying Landau gauge; the 
associated gauge parameter would be the inverse of the ghost form factor at zero momentum, or 
equivalently in our language, the value of cj(0). Moreover, these different gauge-fixing conditions 



should be related [ ]31| , |32| , [33[ ] to different Gribov copies. If this would really be a matter of gauge 
fixing, it looks rather suspicious that an a priori general framework of gauge fixing would make no 
sense in J = 2, where any decoupling, in the sense of ^(0) > 0, is excluded. The possibility of a 
d = 2 decoupling/massive solution would nevertheless be expected in the philosophy of ^,33] 
(see again [ p^ for a critical discussion). 



5. From gluon to ghost propagator 

Armed with the current high precision lattice data for ghost and gluon propagators, it remains 
the question if these can be described in a single analytical framework, not only qualitatively, but 
also quantitatively. An interesting approach to this problem is the philosophy of [|^]: given a 
quantitatively fine description of the lattice gluon propagator, does the corresponding solution of 
the ghost Dyson-Schwinger equation produce an equally fine quantitative description of the corre- 
sponding lattice ghost propagator, this without changing any of the input parameters. It was found 



in [ 34 ] this is not a trivial requirement, as depending on the strength of coupling constant, a set 



^Notice that this does not imply anything for S^{0), as we would then simply prove the rather obvious bound 
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of solutions can be found. Of course, in principle the coupling constant in a given renormaliza- 
tion scheme is a fixed quantity. A similar observation has been made in [||]. Also the role of the 
ghost-gluon vertex cannot be underestimated for quantitative results [|27|]. 

We follow a slightly different route, based on the papers ||3|, ^]. In [ |l0[ ^ ], high quality fits 
were presented to the lattice gluon propagator md = 2,3,4. In the particular case of d = 3,4, the 
used fitting functions are in a 1-1 correspondence with the tree-level propagators of the Refined 
Gribov-Zwanziger (RGZ) formalism [|l^, [l^, thereby fixing all the input (massive) parameters of 
the RGZ approach using lattice input. These masses correspond to condensates, which existence 



is favoured based on dynamical reasons [13|. Said otherwise, they correspond to another vacuum 
state with lower energy. The natural question arising then is whether perturbation theory around this 
nonperturbative vacuum is sufficient to capture all relevant dynamics for e.g. the ghost propagator. 



Work in this direction is currently being undertaken in [35]. Preliminary results do point towards 
the need of including gluon-ghost vertex corrections to match the one-loop ghost propagator with 
its lattice counterpart over the full momentum range. Nevertheless the correct leading infrared 
behaviour is already recovered with one loop perturbation theory. It is also instructive to note 
here that the coupling constant is not a free parameter in this setting, as it ought to be related to 
the renormalization scale using its (one-loop) renormalization group equations in e.g. the MOM 



scheme |35] 
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